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Abstract 

Gas clouds under the influence of gravitation in thermodynamic equi- 
librium cannot be isothermal due to the Dufour efi'ect, the energy flux 
induced by density gradients. In galaxy clusters this effect may be re- 
sponsible for most of the "cooling flows" instead of radiative cooling of 
the gas. 

Recent observations of galaxy clusters with high spatial resolution have 
shown that in most of them the intracluster plasma is far from being isothermal, 
exhibiting strong decrease of the plasma temperature in the vicinity of individ- 
ual large galaxies. This effect is commonly attributed to "coohng flows", that 
means, increasing radiative energy loss with increasing density of the gas, which 
is gravitationally attracted by the galaxies. But several discussions (Voigt et al. 

Markevitch et al. 0) have shown that most of the temperature gradients 
would be washed out by thermal conduction, if the thermal conductivity of the 
plasma were of the magnitude given in the classical book by Spitzer [H]. 

It has been proposed that magnetic fields may reduce the thermal conduc- 
tivity. The existence of magnetic fields in the order of 0.1 — 1/iG has been 
confirmed by various observations (see e.g. Dolag et al. But to produce 

the required reduction, ordered magnetic fields of a still higher magnitude would 
be necessary, as the heat flux reduction works only perpendicular to the mag- 
netic field lines. Narayan and Medvedev have shown that fields fluctuating 
on a large range of length scales can produce only minor changes in the thermal 
conductivity. 

In all these discussions it is implicitly assumed that in thermodynamic equi- 
librium conduction leads to a uniform temperature, though it is well known 
that density gradients may cause a flux of energy and thus induce tempera- 
ture gradients, an effect known from the textbooks as "Dufour effect" (see e.g. 
Hirschfelder et al. 

Speaking of equilibrium means that the net fluxes of all properties determin- 
ing the state of a system are locally balanced. Normally in a gas we associate 
this with constant pressure, density and temperature. But in a system, which is 
influenced by volume forces such as gravitation, the equilibrium state may well 
exhibit gradients of the state variables. Mass flow equilibrium in these systems 
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is obtained by the balance between the gravitational force and the counteracting 
pressure gradient, which is related to the density gradient by some equation of 
state. 

But if there is a density gradient, there must be also a temperature gradient, 
balancing the Dufour effect, to obtain zero net energy flux. Below we will derive 
the equilibrium conditions for a system, which is fully relaxed with respect to 
flux of mass, momentum and energy, within the scope of a simplified kinetic 
model. 

For this purpose we determine the net flux of mass and energy through 
a control area AF under the condition that the particles in the gas have a 
Maxwellian distribution everywhere, but with number density and mean energy 
varying perpendicular to the control area. The net flux is calculated then under 
the assumption, that any particle crossing the control area retains its velocity 
and direction, which it has obtained in the last collision, one mean free path 
from the surface, and in addition is subjected to some acceleration from the 
force field. In the energy balance the gain and loss of kinetic energy due to 
the changing gravitational potential will be omitted, as it cancels out from the 
balance exactly, when the net particle flux through the control area is zero. 

Denoting the direction perpendicular to the surface as z and the angle be- 
tween this and the direction of particle motion as 6, the flux of some property 
G through the surface is 

J AF cos 9{j+ -r)dn, (1) 

where and arc the normal components of the flux of G carried by particles 
moving towards the surface in positive and negative z direction under the angle 
9. The integration has to be taken over half the solid angle 7r/2 > 6* > 0. 
Denoting the amount of property G carried by particles moving with velocity v 
by g{v), the flux components of g{v) are 

j+ ^ig{v)-^^Xcos9)ivcosd + b-) (2) 
dz V 

J- = (g(v) + ^^Xcos9){vcos9-b^) (3) 

A is the mean free path, which may depend also on the velocity. The last term 
denotes the acceleration of the particles during their flight from the last collision 
to the control surface. Integrating eq.(^) over d^l yields the equation 

A ^ dg(v) 2bX , \ 
J = -ttAF f Ai;-^ - —9{v)j (4) 

To determine mass and energy transport by particles moving with velocity 
V, g has to be set to 

gu n{z)mf{v,z)dv, g^ ^ n{z)m—J{v,z)dv (5) 

In case of a neutral gas n{z) and m are the density and mass of atoms, /(tj, z) 
is the distribution function 

f{v)dv^^^e-'^dv (6) 
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P{z) being the abbreviation f3{z) = yj2kT{z) /m. To calculate the quantity 
dg/dz we also need the derivative 



df _d(3 4 /2w4 3w2\ _^_ld/3/2w 



-3)/M (7) 



dz dz ^ J P dz\ P 

While in a neutral gas the parameter b in eq.Q) is the normal gravitational 
acceleration, in a plasma the action of the force is somewhat indirect. Energy 
transport is mediated preferably by electrons. Thus in eqs.© n{z) and m are 
the electron density and mass. But in this case acceleration by the external force 
is acting also on the ions, which transfer the force to the electrons by Coulomb 
interaction. Thus the parameter b in eqs.Q and Q is not the gravitational 
acceleration of free moving electrons. But for the equilibrium condition of zero 
transport the absolute magnitude of b is not relevant, when we only want to 
know, under which conditions the net flux of electrons and the flux of energy 
carried by them are in balance. 

An additional difference between neutral gas and plasma results from the fact 
that the collision cross section between neutral atoms is nearly independent of 
the collision energy, so that the mean free path A does not depend on v. In 
plasmas, where Coulomb interaction is dominant, the mean free path increases 
with . Thus in the further calculations we set A ~ Ao(w/wo)" with a = for 
neutral gas and a = 2 for a fully ionised plasma. 

Introducing eqs.© to (TJ into eq.Q), we get for the mass and energy flux 

iM = 2/3" ( xp^ + x(2x^ -3)n^~^-] h{x)dx (8) 
\ dz dz p X J 

jE = /3"+2fx3/3^+a;3(2a;2-3)n^-^x^ h{x)dx (9) 
\ dz dz P J 

with x^v/P and h{x) = 2^mAF{Xo/vg)x'^+°'e^''\ 

Integrating these equations over all x, we obtain the conditions for zero mass 
and energy flux: 

P^F3+^ + [2F5+„ - 3F3+^] - = (10) 

dz dz p 

f^^F^+c + [2F7+„ - 3F5+a] - ^Fs+a = (11) 

The abbreviation F^ stands for the integral Fj, — x^e~^^ dx. The differ- 
ent arguments of F^ in eas. HlUI) and Hll|l can be eliminated by the recurrence 
formula Fk+2 = (fc + l)/2 x Fk, so that we finally obtain the equilibrium condi- 
tions for mass and energy transport: 

2 + a^dn (2 + a)(H-a) dP 2bn ^ 
^^d^ + 2 " - — ^ ^ 

4 + a^dn (4 + a)(3 + Q;) dP 2bn ^ 

— 7^ — P—, — I TT-—, 77- — (lo) 

2 dz 2 dz P ^ ^ 
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It is immediately evident that this is a set of hnear equations for the gradients 
dfi/dz and dn/dz. If there are no external forces (5 = 0), there exists only the 
trivial solution dp/dz — and dn/dz — 0. 

In the presence of volume forces such as gravitation eas. (|12|l and are in- 
homogeneous and allow a non-trivial solution, which constitutes a fixed relation 
between temperature gradient and particle density gradient or, because of quasi- 
neutrality, also between temperature and mass density gradient. Eliminating b 
leads to the relation 

dn , , dB , , 

/3— + 2a + 5 = (14) 
dz dz 

with the solution 7i/3^^"+^^ — const. Replacing the number density by the mass 
density p and (3 by ^2kT jm we finally find the condition pT'°l^ — const, for 
neutral gas and pT^/^ — const, for a fully ionised plasma. 

Though the relative temperature gradients are small compared to the asso- 
ciated density gradients, they may well be important in the hot corona of stars 
or in the intergalactic gas in galaxy clusters, where temperature changes by a 
factor of three are reported (Voigt et al. [7|), while the gas density may change 
by a few orders of magnitude in the vicinity of large galaxies. It may even 
be that most of the observed temperature drop in the "cooling flows" around 
these galaxies is not caused by radiative cooling but by the " conductive cooling" 
associated with the density gradients. 

It should be noted that due to the Dufour effect complete galaxy clusters 
would have cooled off, if thermal conductivity were at the Spitzer value, unless 
the clusters are embedded in a very hot but thin intergalactic plasma. The 
existence of a diffuse high energy x-ray background (Boldt P) is a strong hint 
to its existence. Matter expelled in the jets of quasars may be the origin of this 
plasma. 
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